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Abstract 

Motivated by recent experimental results in ultra-thin helimagnetic Holmium films, we have 
performed an extensive classical Monte Carlo simulation of films of different thickness, assuming a 
Hamiltonian with six inter-layer exchange constants. Both magnetic structure and critical proper- 
ties have been analyzed. For n > 16 (n being the number of spin layers in the film) a correct bulk 
limit is reached, while for lower n the film properties are clearly affected by the strong competition 
among the helical pitch and the surface effects, which involve the majority of the spin layers: In 
the thickness range n = 9 -f- 16 three different magnetic phases emerge, with the high-temperature, 
disordered, paramagnetic phase and the low-temperature, long-range ordered one separated by 
an intriguing intermediate-temperature block phase, where outer ordered layers coexist with some 
inner, disordered ones. The phase transition of these inner layers displays the signatures of a 
Kosterlitz-Thouless one. Finally, for n < 7 the film collapse once and for all to a quasi-collinear 
order. A comparison of our Monte Carlo simulation outcomes with available experimental data is 
also proposed, and further experimental investigations are suggested. 

PACS numbers: 64.60.an,64.60.De,75.10.Hk,75.40.Cx,75.70.Ak. 
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I. INTRODUCTION 



Surface and nanometrical objects are important both for their possible implementation 
in the current technology and from basic research point of view. In collinear magnetic thin 
films intriguing behavior, as transition temperature depending from the film thickness — n— 
or critical exponent crossover, were observed 1, 2, ^, Q, S, s| • 

Nowadays, the most fervent interest has moved towards frustrated systems, where a non- 
collinear order, characterized by a possibly large modulation, is established. Some rare-earth 
elements (as Holmium, Dysprosium or Terbium) and their compounds are typical examples 
that the nature makes available to investigate such peculiar behaviours, in view of the variety 
of magnetic arrangements, as helix, spiral or longitudinal-wave, that can be observed in bulk 
samples of such materials. M. Further examples of helicoidal structures can also be met in 
multiferroic materials . 1 1 Of ] and itinerant systems, as MnSi and FeGe [121 ]. 

In magnetic systems with frustration, the lack of translational invariance due to the 
presence of surfaces can result especially important for ultra-thin film samples where the 
thickness is comparable, or even lower, with the wave length of the ordered magnetic struc- 
ture observed in the bulk. It is worthwhile observing that when these conditions are met 
a sweeping change of the magnetic structure behaviour could be found. Many fundamental 
features related to such systems have not yet been exhaustively investigated and completely 
understood, and ultra-thin films of rare-earth elements are still among the most intriguing 
layered systems to be studiedflja]. 

From an experimental point of view, the availability of new sophisticated growth and char- 



acterization techniques 



14J ] has allowed to extensively investigate the properties of such mag- 



netic nano-structures. For instance, interesting experimental data on thin films of Holmium 
(whose bulk samples show helical order along the c axis, perpendicular to film basal planes) 
were obtained 13, UM, ll6|, LlTJ by neutron diffraction and resonant soft X-ray experiments. 
By looking at the static structure factor S(Q) (Q = (0, 0, Q z ) being the wave vector of 
the incommensurate magnetic modulation along the film growth direction z), it has been 
shown that the critical behaviour of Holmium thin films markedly differ from that of films of 
transition metals, which usually display a collinear ferroma gnetic ( FM) or antiferromagnetic 



13 



15 



16 



17 identified a thickness 



(AFM) order in the bulk. Mainly, the authors of Refs. 
uq ~ 10 mono-layers [ML] (comparable with the helix pitch of bulk Holmium ~ 12 ML) 
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which was interpreted as a lower bound for the presence of the helical ordered phase. In 
fact, they observed that the transition temperature dependence on film thickness does not 
follow the usual asymptotic power law [3], but rather an empirical relation 

T N (oo) - T N (n) 1 , 

T N (n) ~ (n - n Q y ' 1 j 

can be devised [3], where T N (oo) and T N (n) are the ordering temperature of the bulk system 
and of a film with thickness n, respectively, while the exponent A' has not an universal 
character. Interesting enough, it appears that this empirical relation is not peculiar of 
helical-like structures only but it results more general, being observed in other ultra-thin 
structures characterized by a magnetic modulation as well: An important example is given 
by Chromium films, where at low temperatures an incommensurate spin density wave is 



present 



18|- 



13 



and 



15 in order to un- 



A mean field approximation (MFA) was proposed in Refs. 
derstand the experimental outcomes from Holmium films: MFA allowed to obtain a first, 
rough estimate of the threshold thickness n defined in the empirical relation (pQ), but it 
also revealed that for thicknesses close enough to n Q the paramagnetic and helical phase 
can be accompanied by a more complex block phase, where groups of ferromagnetically 
ordered layers pile-up in an alternating antiferromagnetic arrangement along the c axis. 
As it is well known, the MFA completely neglects thermal fluctuations, which are however 
strongly expected to play a fundamental role in the critical behaviour of low-dimensional 
magnetic systems: therefore, adding thermal fluctuations not only give strong quantitative 
adjustments of the MFA estimates of critical quantities as T^{n) or n , but could also make 
unstable some ordered structure as the block phase met above. 

In order to overcome such issues and deepen our understanding of critical phenomena in 

Holmium ultra-thin films, we performed extensive classical Monte Carlo Simulations (MCS): 

□ 

Preliminary results already showed [19( that thermal fluctuations do not destroy the block 
phase, which instead acquire a much richer structure, with disordered inner layers interca- 
lating ordered ones and undergoing a Kosterlitz-Thouless (KT) phase transition 36J as the 
temperature lowers. 

A complete account is here given of the results of our simulations for film thickness in the 
range n = 6 . . . 36 ML. The paper is organized as follows: In Sec. [IT] we shall briefly recall 
relevant properties of Holmium, and introduce the magnetic model Hamiltonian. Sec. II III 
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is devoted to describe the Monte Carlo method and the estimators employed to evaluate 
the physical quantities relevant for magnetic films with non-collinear order. In Sec. IIVI the 
Monte Carlo results about the magnetic order establishing at low temperature are shown for 
different thicknesses. The role of thermal fluctuations is discussed in Sec. |Vj In particular, 
Sec. IV Al will report a detailed study of the temperature regions where the single layers 
display a critical behaviour, the structure factor close to these regions being deeply analyzed 
too, given its fundamental relevance in an experimental mindset; Sec. IV Bl will be devoted to 
the global film properties. All the results reported in the previous sections will be compared 
and discussed in an unifying framework in Sec. IVIl where we shall also gather our conclusions. 



II. MODEL HAMILTON! AN 



The magnetic properties of Holmium have been intensively investigated both experimen- 
tally and theoretically [2]. The bulk crystal structures is known to be hexagonal close- 
packed (hep). The indirect exchange among the localized 41 electrons manifests as an 
RKKY long-range interaction of atomic magnetic moments; the experimental data about 
the low-temperature magnetic moments arrangement in Ho can be reproduced assuminga 
FM interaction between nearest neighbor spins lying on the ab crystallographic planes 3] , 
while along the c crystallographic axis interactions up to the sixth neighboring layers must 



be allowed for (see for example Ref. 120). It's just the competing nature of the latter that 
below Tjy(oo) = 132 K gives rise to an incommensurate magnetic periodic structure, which 
can be modeled as an helical arrangement of the magnetic moment vectors along the direc- 
tion (henceforth denoted as z) parallel to c, i.e. perpendicular to the ab crystallographic 
planes, where the magnetic vectors prefer to lye as a consequence of a single-ion easy-plane 
anisotropy. The average local spin vector at low temperature can thus be expressed as: 

Si = s(fi) = (s± cos \Q- ri),s± sin (q • fy ,s z } , (2) 

where Q = (0, 0, Q^ lk ) is the helical pitch vector, and cQ b z ulk ~ 7r/3, c being the lattice 
constant |7(] along the z direction (i.e., c/2 is the distance between nearest neighbouring ab 
spin layers). In addition, the crystal field bring into play other different kinds of anisotropies 
that, at temperatures well below T/v(oo), are able to change the helical shapes in conical 
ones, or force the magnetic structures in a bunched helix which is commensurate with the 
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lattice [2l|,l22|. 

Theoretical investigations have shown that incommensurate magnetic bulk structures 
(observed, besides Holmium, also in other rare-earth elements, as Dysprosium and Terbium) 
can be well obtained by a MFA |23j from a simple Heisenberg model with only three coupling 
constant, the first one, Jo > 0, describing the FM in-plane interactions, while J\ and J 2 are 
the effective coupling between ions on neighbouring (NN) and next-neighbouring (NNN) 
planes, respectively. Whatever the sign of J±, the MFA finds a helical structure when 
J 2 < 0, i.e. AFM, and the conditions | J 2 | > | Ji|/4 is met. 

It is worthwhile to recall that when dealing with ultra-thin films the assumption of being 
allowed to retain the same Hamiltonian able to describe the bulk structure is absolutely not 
guaranteed to be correct: Indeed, real film samples can be strongly affected by defects, strain, 
thickness uncertainty (~2 ML) or interaction with the substrate 24( (typically Y/Nb or 
W(110)). The latter can be particularly relevant, as it can change, sometimes dramatically, 
the single ion anisotropy and the strength of the interaction constants with respect to bulk 



samples. Furthermore, the lack of inversion symmetry can bring into play t 
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re antisymmetric 



and 



271 for perovskite 



Dzyaloshinskii-Moriya (DM) interaction [25[ (see for instance Ref. 
multiferroic RMn03 with R = Gd, Tb, or Dy) and possible surface anisotropics for Dy/Y 
multilayer films [281 ]. 

While always remembering such possible drawbacks, in our investigation of Ho thin films 
we employ a Heisenberg model Hamiltonian which has proven useful to describe Holmium 
bulk samples. We thus define: 



N 



N 



n 



1,3 



k=l 



(3) 



L being the lateral 



where N = the total number of magnetic ions, L r = L 

film dimensions and n the number of layers (see Fig. [Th), and <7j are classical unitary vectors 
representing the total angular momentum of Ho ions (i.e. a, = Ji/\Ji\). As a consequence, 
the magnitude of the angular momentum of Ho ions is embodied in our definition of the 
coupling constants appearing in Eq. ([3]) as = |] l 7'|(| l 7'| + 1) J^-, with |j7i| = 8. For the 



20, while the 



inter-layer exchange parameters J,;,- we have assumed the values given in Ref. 
easy plane anisotropy D z = 16 K [23]. 

In Fig. [Tb we show a schematic representation of the hep lattice structure and of the 
exchange interactions included in the model Hamiltonian ([3]) . On the hexagonal basal planes 



(a) 



helix 
z-direction 



t 



n 




(b) 




FIG. 1: (color online) (a) Film geometry: N = n x L x x L y is the total number of spins, n 
is the film thickness, i.e. the number of spin layer (free boundary conditions are taken along z- 
direction), and L x ^ y are the layer dimensions (periodic boundary conditions are applied along x and 
y directions), (a) Schematic representation of hep Ho structure and exchange interactions: only Jo 
(in-plane, green lines), J\ (NN planes, red lines), J2 (NNN planes, blue line), and J3 (cyan lines) 
are represented. The numerical values of the coupling constants employed in our simulations are: 
J = 55.68 K, Ji = 24.9168 K, J 2 = 5.4288 K, J 3 = -3.3686 K, J 4 = -4.5936 K, J 5 = -0.2784 K, 
J fi = -2.5056 K. 



only a NN, FM interaction Jo > (green lines) is considered; along the c-axis we instead 
allow for interactions up to the sixth neighboring layer, with a total coordination number 
C = 30. 

Sign and magnitude of the out-of-plane interaction constants Ji ... J 6 were determined 



in Ref. |20j in order to reproduce, with the correct pitch vector, the helical ground state 
along the c-axis observed in the experimental data. However, at the best of our knowledge, 
also neutron scattering experiments investigated the dynamical properties of Holmium only 
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along the c-axis 22j, so that a direct measure of the in plane FM coupling constant is still 
lacking, and only mean field estimates are available, which set it at about 300/zeV [l^. This 
allows us to consider J as an almost free fit parameter to be adjusted in order to fix the 
correct value of experimentally accessible quantities. By Monte Carlo simulations we find 
the value T^ cs (oo) ~ 124 K for the bulk ordering temperature by setting J ~ 400 /zeV: 
bearing in mind the cautions given above about the possible quantitative difference among 
bulk and film samples, this is the value of Jo we have used in all the following simulations of 
thin films, without attempting any further possibly meaningless quantitative adjustment. 

III. MONTE CARLO METHODS AND THERMODYNAMIC OBSERVABLES 

Our study of the magnetic properties of thin rare earth films was done by extensive 
classical MCS. Thickness n from 6 to 36 and lateral dimensions L x = L y = L = 8 . . . 80 
have been analyzed. As we are working on film structures, free boundary conditions in the 
thickness direction z are obviously taken, while the usual periodic boundary conditions are 
applied in the L x L planes (see Fig. Wfr), which coincide both with the ab crystallographic 
planes and the easy plane for the magnetization. 

Simulations were done at different temperatures; the thermodynamic equilibrium is 



reached by the usual Metropolis algorithm 



29( and over-relaxed technique 30]. The latter 



was employed in order to speed-up the sampling of the whole spin configuration space: In- 
deed, the competitive nature of the exchange interactions in the Hamiltonian ([3]) and the high 
coordination number lead to a long time needed to reach thermal equilibrium, especially in 
the critical region. We have thus resorted to a judicious mix of Metropolis and over-relaxed 
moves in order to reach the goal in a reasonable time. Usually, one "Monte Carlo step" is 
composed by one Metropolis and four/five over-relaxed moves per particle, discarding up to 
5xl0 4 Monte Carlo steps for thermal equilibration; at least three independent simulations 
where done for for each temperature. 



The Monte Car 
togram technique 



o data analysis have benefited from the employment of the multiple his- 



31 



321 ]. which allows us to estimate the physical observables of interest 
over a whole temperature range by interpolating the results obtained from simulations per- 
formed at some chosen, different temperatures. The outcome of the method is an estimate 
of the density of state p(£ ) at energy £ obtained by weighting the contributes Pi(£) due to 
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independent simulations made at inverse temperature $ = l/fc^Tj). The independent 
simulations have to be sufficiently close in temperature, i.e. the temperature step must be 



Z„ = X>(£)e-* = £ 



(4) 



chosen roughly proportional to the square root of the inverse of the heat capacity 31] • We 
thus have computed the partition function Zp at any (3 in the range of interest by solving 
iteratively the equation 

-l 

C i,s L fc 

where = f3 — j3k and is the number of independent samples of energy for the k- 
simulation. The index % refers again to the single simulation, while s denotes the energy 
sampling intervals in the i-th simulation. 

As mentioned above, the high number of exchange interactions makes difficult the es- 
timate of the density of state, especially close to the critical temperature, and variables 
like specific heat or susceptibility are extremely sensitive for these systems. Anyway, this 
obstacle has been successfully overcome making use of the estimator, 

c v = ^({e-(S) p ) 2 ) p (5) 
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Iterating several times the 



in the histogram reweighting techniques, as suggested in Ref J 
multiple-histogram algorithm, we have also obtained the variance of the interpolated data 
by bootstrap resampled method, by picking out randomly a sizable number of independent 
measurements mk (between 1 and 5xl0 4 ), and iterating the re-sampling at least one hundred 



times 



31 



3J. 

As we are interested in the phase transitions of Holmium films, it is worthwhile to observe 
that the study of films described by the Hamiltonian ([3]) entails a wide number of funda- 
mental issues. First of all we must consider: i) the intrinsically two-dimensional (2d) nature 
of such magnetic structures; %%) the presence of different interactions, which turn out to be 
FM on the layers (with a SO (2) symmetry) and competitive along the thickness direction n, 
with a possible helimagnetic (HM) order at low temperature (i.e. a Z2 xSO(2) symmetry 



351]); and Hi) the implementation of different boundary conditions for n and L, respectively, 
above introduced. 

About the first issue, it is well known that the critical behaviour of an ideal easy-plane 
magnetic film with continuous symmetry and short-range FM (or AFM) interactions per- 
tains to the 2d XY universality class, displaying a KT behaviour [36] at a finite critical 

8 



temperature. In particular, a crossover from 3d to 2d behaviour is expected when the corre- 
lation length saturates the film thickness. However, from MCS point of view, it may result 
quite difficult to realize such conditions. Indeed, even for large, but still finite, L values a 
sharp transition can not be observed, making it possible to define a 3d pseudo critical point, 
as extensively discussed by Janke and co-workers in Refs. 



37 



38 



Turning to the second issue, in a quasi-2d magnetic system with a continuous symmetry 
the introduction of competing interactions along the direction perpendicular to the film slab 
brings into play the presence of two (in-plane and out-of-plane) correlation lengths, with 



19, under these 



a rather dissimilar behaviour in the critical regime. As analyzed in Ref. 
conditions some new and interesting critical phenomena can be observed. Systems with 
discrete symmetries which present two different correlation lengths were already discussed 
in literature, see, e.g., Refs. |39| and Q. 

Moving to the last issue, we must first of all observe that the identification of a suitable 
order parameter to study the critical properties of non-collinear thin films requires a careful 
analysis of some of their peculiar features. A first trouble is the intrinsic difficulty represented 
by a helical order parameter associated to a wave vector Q: In the bulk system the virtually 
infinite size of the system, summing up an infinite number of in-phase contributions, makes 
a clear peak emerge at wave-vector Q in the static structure factor in the ordered phase. 
In films, the presence of broad peaks is on the contrary expected in a wide temperature 
range as, a consequence of the intrinsic finite size nature of the system, thus jeopardizing 
the identification of a well defined critical temperature Tjv(n) from the sole analysis of the 
peaks appearing in the structure factor. Secondly, as we will discuss in the next Sections, 
a naive structure factor analysis could be not enough to distinguish between an HM order 
and other ordered phases that can be present. For these reasons, it is necessary to resort 
io other observables related to the HM order. A first choice can be found in the chirality 
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421 ]. which can be defined on film as: 

1 



~~ 3(n-l)L 2 smQ b z ulk - 

z li 



X °l+l,iY > ( 6 ) 



where I labels the planes, starting from one of the two film surfaces, and i locates the spin 
on the plane. As we shall see in the next Sections k represents quite a good quantity to 
locate the critical temperature for the HM phase. In view of point i) discussed above, it is 
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useful to introduce a FM order parameter for each layer /: 

M t = yj{MfY + (Mf) 2 + (M? ) 2 , (7) 
where Af£_= A E; ^ with 

a = x } y } z, and consequently the average order parameter of 

the film [42, 

1 n 



M 



(8) 



It is worth to observe that the order signalled by the quantities defined in Eq. (J7J) and (JHJ) 
do not directly entail the existence of an HM or fan structure in the film. 

The critical nature of the observables defined in Eqs. (jSJ), 0, an d dU), is better revealed 
by looking at the following derived quantities (O = k,M[, M): 



(Xo) = N(3 ((O 2 ) - (O) 2 ) , 



(9) 
(10) 

(11) 



which, at the critical temperature, display a peak that can be characterized by the usual 
finite size scaling theory 45(. In particular, for large enough L, Tjj(n) approximately scales 
as 

T^H^T-H + CL- 1 ^, (12) 

where C is system dependent constant, while v is the correlation length critical exponent. 
A final quantity we employ in our investigation of the critical properties of films is the 



Binder cumulant 



46 



47|: 



Ua 



(Q 4 ) 



(13) 



which allows us to locate the critical temperature by looking at the intersection of the graphs 
of U as a function of T obtained at different L: At T = Tjy(n) such crossing becomes a 
"nontrivial fixed-point" 46|, |47|]. Moreover, one can examine the ratios Uli/Ul (for sizes L 
and V > L) as temperature function, looking for a unitary ratio at the critical temperature. 



IV. MAGNETIC STRUCTURES AT LOW TEMPERATURE 



In this section we will present and analyze our Monte Carlo results for the overall magnetic 
behaviour of film samples of different thicknesses, from a bulk-like structure with n = 36 to 

10 




FIG. 2: Normalized magnetic vector profiles of each layer at T = 10 K, layer dimensions L X)V = 80 
and different film thickness n: (a) n = 36; (b) n = 16; (c) n = 12; (d) n = 11; (e) n = 10; 
(f) n = 9; (g) n = 8; (h) n = 7; (i) n = 6; the error bars are included in point dimensions. 
Inset in (a): the magnified zone of the main graph (here shown for the 6 th , 18 th , and 30 th planes 
from left to right along mf, respectively) clearly shows that the helix pitch and lattice structure 
are incommensurate; inset in (b): the magnification shows the magnetization of the first surface 
layers: 1 st , 14 th , 2 nd , 15 th , 3 rd , and 16 th from top to bottom along the m\ , respectively. 

a very thin film of n = 6 layers, at temperature T = 10 K, i.e. well below T N (n); the lateral 
dimension of the films is taken constant at L = 80, having checked that at this temperature, 
far from the critical region, such value of L well represents the thermodynamic limit for all 
practical purpose. 
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FIG. 3: (Colors online) Angle A(f>i = cfti+i — <pi between magnetic moments in NN layers (I + 1, 1) 
at T = 10 K, and L = 80, for representative thicknesses: n = 16 (black circle), n = 12 (green 
diamond), n = 9 (red square), and n = 6 (blue up-triangle) . Error bars included in the point 
dimensions. 

In Fig. [2] the normalized in-plane magnetic vector (jnf,m^) profile of each layer I is 
reported: For thicknesses greater than 12 ML (which coincide roughly with the helix pitch 
of bulk Holmium) a behavior essentially unaffected by surface effects is observed in almost 
the whole sample, with a typical HM order (Fig. (2^-b, i.e. n = 36 — 16, respectively). 

As emphasized in the the inset of Fig. [2b the magnetizations of NN planes close to the 
surfaces form angles well lower than those observed in the bulk, an expected consequence 
of the increasing lack of interactions on one side of the planes as the surface is approached. 
Such effect can be better looked at by defining the magnetization rotation angle between 
NN planes A<pi = <pi + i — <pi. In Fig. [3] A<pi is displayed for some representative values of 
the thickness: for thick samples, surface effects are especially strong only on the first three 
layers on each film side, and this explains why while for n > 9 an almost bulk behaviour can 
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FIG. 4: (Colors online) Structure factor S(0,0,Q z ) (continuous line) vs. Q z at T = 10 K and layer 
dimensions L x ,y = 80 for different film thicknesses: (a) n = 36; (b) n = 16; (c) n = 12; (d) 
n = 11; (e) n = 10; (f) n = 9; (g) n = 8; (h) n = 7; (i) n = 6. Red dashed lines and green 
dot-dashed lines are the structure factor components along the x or y spin space directions, see 
text. Q z is measured in reciprocal lattice units 2ir/c 



be observed, at least for some inner planes, the scenario changes significantly when n drops 
below 9 (Fig.[2j-h). 

The characterization of the magnetic order can be further pursued by looking at the 
static structure factor S(Q) (where Q = (0, 0, Q z )), i.e. to the Fourier transform of the 
spin correlation function along the z-direction of the films. S(Q) is reported in Fig. H] 
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(continuous line), together with its in-plane components S XX (Q) and S yv {Q) (red and green 
line, respectively; x and y directions do not obviously have any special meaning and can 
be chosen at will: here we use the same orientation already employed in Fig. [2]). Once 
again, for n > 11 (Figs. Hh.-d) both the global structure factor and its components shows a 
clear peak at Q™ ax ~ 0.17, a value in total agreement with the bulk one Q b z ulk ~ 1/6. On 
the contrary, for n ^ 10 a fan-like structure appears, signalled by the emergence of an FM 
component (i.e. a maximum at Q z = 0) in S XX (Q) or S yy (Q) (Fig. Hf-h), while for n < 7 
a quasi-collinear spin arrangement is finally reached, as testified by the single maximum of 
S(Q) itself at Q™ ax = (Fig. H). 

The results discussed so far show that, at low temperature, the progressive film thickness 
reduction does not seem to lead to a sudden helical order suppression, but rather to induce 
a gradual passage to a fan-like order associated with a helix distortion due to the surface 
effects, until a permanent collapse to an almost collinear order occurs for n < 7. Summing 
up we can roughly assume that MCS data show that for thickness n < 9 the helical order is 
substantially absent. 

Such results can be considered in fairly agreement with the experimental outcomes: In 



fact, in Ref. 



15l the authors identified the thickness no — 10 ML as the value indicating the 



complete lack of helical order, and a thickness uncertainty around about 2 ML 24j must be 
taken into account. 



V. MAGNETIC STRUCTURES IN THE ORDER-DISORDER BOUNDARY RE- 
GION 

A. Layer's Magnetic Behaviour 

This subsection is devoted to the investigation of single layer critical properties for some 
n values. At the beginning, our attention will be focused on thicknesses close to the Ho helix 
pitch (12 ML). For this purpose, the order parameter for each spin layer, as defined in Eq. (JTj), 
is evaluated, together with its Binder cumulant, Eq. (|T3j) . and its susceptibility, Eq. (Q. 
Hereafter, we will denote with the symbol Tc n (/) the transition temperature of the I— th 
layer of the film of thickness n. 

Susceptibility and Binder cumulant for the first six layers are plotted in Fig. as a 
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FIG. 5: (color online) Susceptibility (a), and Binder cumulant (b) vs. temperature, of the order 
parameter Mi for thickness n = 12 and L = 12 (black circle), 16 (red square), 24 (green diamond), 
32 (blue up-triangle) , 40 (purple down-triangle), and 48 (brown star). I runs from the first to the 
sixth spin layer. Close to the critical temperature To^ijt), the continuum lines are obtained by 
multiple histogram technique. 

function of temperature for different values of lateral dimension L at n — 12. A critical 
region in a wide temperature range around T ~ 120 K is observed in Fig. [5^ for all planes 
but the central ones (the 6 th and, for symmetry reason, the 7 th ) which are definitely still in a 
paramagnetic state, displaying instead a critical region shifted around T ~ 114 K. Using the 
Binder cumulant, Eq. ([To 7 ]) , for different values of L we can estimate the single layer transition 
temperature ^12(6) = 113.4(4) K of the inner planes and Tc7 5 i 2 (l . . . 5) = 120.3(4) K of the 
external ones. 

The intriguing landscape here observed for n = 12 is present in the whole range 
9 < n < 16. A summarizing picture of the single layer critical temperature Tc7 n (/) vs. 
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FIG. 6: T c , n (l) vs. layer index I for n = 20 (a), 16 (b), 12 (c), 9(d), and 8(e). The vertical dashed 
lines denote the position of the bisecting plane of the film, beyond which single-layer properties 
repeat by symmetry. 

plane index I for n = 20, 16, 12, 9 and 8 is given in Fig. [6) For the thicker film here an- 
alyzed (n = 20, Fig. [6^) 7c-2o(0 * s the same f° r every layer, and coincides with the es- 
tablishing of HM order in the film, as expected for the bulk system, where the critical 
temperature can be obtained both through the chirality, Eq. ([6]), and by Eq. (jSj). For 
n — 16 (Fig. [6b) we observe a structure more complex than that we find in the films 
with n = 12 and n = 9 (Figs. [6b and Fig. Eli respectively). Indeed, as discussed in 
the 5 th , and 6 th (and, by symmetry, the 11 th and 12 th ) planes lose 



details in Ref. 
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their order at a lower temperature, 7b jl6 (5, 6, 11, 12) = 117.8(2) K, than the others, where 
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FIG. 7: (color online) Angle formed by the magnetization vectors of NN planes, for n = 12 and 
L = 48 at different temperature: T = 10 K (solid black circle), T = 110 K (open green diamond), 
and T = 118 K (open red square) [48]. Error bars lye within point size. 

T Cil6 (1...4,7...10,13...16) = 120.1(4) K. 

In order to better understand the magnetic structure in these temperature ranges, we 
examine the T-dependence of the magnetization rotation angle A<fri. For the sake of clarity 
in Fig. [7] the system n = 12 is again analyzed. When Tc 5 i 2 (6) < T < Tc t u(l ... 5) (see, e.g., 
red square and line in Fig. [7j) the ordered layers distinctly display a block structure where 
A(pi among the first (last) five planes is almost zero, i.e. < 10°: at the same time the angle 
formed by the magnetization of the two blocks is about 180°j48|. Only for T < T<7i 2 (6) 
(black and green symbols and lines in Fig. [7]), the function A<pi displays the expected thin 
films helimagnetic behaviour discussed in Sec JIVI (for comparison, see Fig. 5a of Ref. [lol 
where the same quantity at T = is discussed within MFA). 

We can graphically represent the block magnetization arrangement for n = 12 as f i | i | 
00 II!!!, where the circle represents the disordered planes, and the arrows the ordered ones. 
As above anticipated, the spin block phase is obtained down to n = 9 (Fig. [6jl), where we 
get an arrangement tttt ° III!) an d up to n = 16 (Fig. where a much more intricate 
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FIG. 8: (color online) Temperature evolution of 5(0,0, Q z ) vs. Q z . (a): thickness n = 12 and 
L = 48 at T = 10 K (black line and circle), 113 K (red line and square), 116 K (blue line and 
triangle), and 121 K (green line and diamond). Dashed black line: S{Q) of the saturated block 
structure TTTTT 00 IIJ.il- Both the black curve (T = 10 K) and the dashed one have been divided 
by a factor equal to ten. (b): thickness n = 8 and L = 48 at T = 10 K (black line and circle), 50 K 
(red line and square), 90 K (blue line and up-triangle) , 105 K (green line and diamond), and 110 K 
(orange line and down-triangle); all temperatures are lower than Tcs(l ... 8). 

layout, i.e. TTTT 00 llll 00 TTTT; is observed. It is worthwhile to observe that the AFM 
alignment of consecutive ordered blocks reveals as the medium range alternating inter-layer 
exchange coupling give rise to an effective AFM interaction between blocks. 

A further insight in these block phases, especially relevant from an experimental point 
of view, is obtained by analyzing the behaviour of the structure factor close to Tc n (l). 
In Fig. [8^ S(Q) for n = 12 in a wide temperature range is plotted. In particular, one 
temperature value just below Tc,i2(6) (T = 113 K, red line and square), one in the block 
phase region T C ^ 2 {Q) < T < T<7 12 (1. . . 5) (T — 116 K, blue line and up-triangle), and one 
just above Tc,i2(l ... 5) (T — 121 K, green line and diamonds) have been chosen. As already 
observed in Sec. [TTJ, the prominent broadening in the whole temperature range of the peaks 
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FIG. 9: (color online) Logarithm of the susceptibilities xm 5 (square) and xm (circle) versus the 
logarithm of the lateral film dimensions L (error bars included in the symbols), for n = 9, at 7^9(5) 
and Tc 5 19(1 ... 4, 6 ... 9), respectively. 

displayed by S(Q) is mainly a consequence of the intrinsic finite size nature of the films. 
Excluding an obvious intensity reduction for increasing temperature, one can immediately 
observe that both shape and peak position are almost unchanged moving from block phases 
(e.g., T = 116 K data) to HM order (e.g., T = 113 K data). A further confirmation of the last 
statement can be achieved by the comparison, again proposed in Fig.[8k., between the Monte 
Carlo outcome for S(Q) at T = 10 K (black line and circle) and the Fourier transform of the 
static block structure T T T T T 00 11111 with saturated magnetization for each ferromagnetic 



plane (black dashed line). As formerly observed for the ultra-thin film n = 16 [191 ]. even for 
n = 12 the two plots have the same peak position and width; moreover they have comparable 
intensities too. We are thus led to conclude for the substantial impossibility to distinguish 
between block phases and helimagnetic order by looking at the structure factor, being it 
able to give information about the global structure modulation only. 



As already observed in Ref. 
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for n = 12, we find that whenever a blocked phase temper- 



19 



ature range occurs, the spins lying on disordered layers are seen to feel a local magnetic field 
due to inter-layer interactions much smaller than that acting on spins on the ordered layers, 
so that they behave as being effectively decoupled from the other ones and display the char- 
acteristic features of a two-dimensional magnet. The different effective dimensionality of the 
critical behaviour of lowest-temperature ordering layers from highest-temperature ordering 
ones, is illustrated in Fig. [9j where an accurate finite-size scaling analysis of the layer mag- 
netic susceptibility xm 5 and of the global susceptibility xm, at Tc,9(5) at Tc i9 (l ... 4, 6 ... 9), 
respectively, is reported for n = 9: making use of the usual scaling relation at the critical 
temperature x ^ L 1 ^, where 7 and v are the critical exponents of the susceptibility and 
correlation length, respectively, the value 7/1/ = 1.72 (4), is obtained from the best fit of xm s 
data, while 7/1/ = 1.94 (3) is the result of the fit of xm- The former value is completely con- 
sistent with the Kosterlitz-Thouless behaviour expected in an isolated two-dimensional, easy- 
plane magnet, while the latter clearly indicates a planar three-dimensional-like trend 371. 138] 
for the system made by the planes 1 -j- 4 and 6 -v- 9. 

We now move to discuss the MC results obtained for n = 8. From Fig. [6fe the lack of 
the ordered-disordered blocks mixed structure at intermediate temperature is apparent, as a 
transition temperature common to all planes (T^^l . . . 8) = 110.6(2) K) is found. In Fig. [8b 
S(Q) vs. Q z for different T is shown: Close to Tc$Q- ■ ■ -8), we have Q™ ax = 0, signalling 
the presence of a FM-like collinear structure. Subsequently, as the temperature lowers, and 
the FM spin arrangement opens towards a more stable fan structure, S(Q) develops a peak 
at Q™ ax = 0.095, but still with a very strong contribution also at Q z = 0. 

The evolution of the structure-factor peak position with temperature is better illustrated 
in Fig.[10j where Q™ ax vs. T is plotted for some significant values of film thickness. For n = 8 
the clear jump from a collinear structure to a fan-like one (reached at T w 90 K) is observed: 
This shows that the onset of order in every plane is by itself not necessarily enough to 
generate the fan structure observed at low temperature. On the contrary, for thickness values 
close to the helical pitch and above (where Q™ ax is essentially independent of temperature) 
the completion of planes ordering, with the transition of inner layers, marks also the onset 
of the overall helical/fan arrangement, while for small n a ferromagnetic alignment again 
stabilizes as soon as the layers simultaneously order. Therefore, the peculiar behaviour of 
Q™ ax {T) for n = 8 can be reasonably attributed to its representing the borderline between 
helical/block ordered structures and substantially ferromagnetic ones. 
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FIG. 10: (color online) Q™ ax vs. temperature for some thicknesses: n = 16 (black line and triangle), 
9 (red line and circle), 8 (green line and square), and 7 (blue line and diamond). For n = 8 the 
stability of the fan phase with respect to the collinear one is reached only below T ~ 90 K (see 
text). 

In view of the previous discussion, we can conclude that the existence of an intermedi- 
ate temperature region, characterized by the presence of spin block structures, in between 
the paramagnetic and the helical ones, seems to be a peculiar feature of non-collinear Ho 
magnetic films with thickness close enough to the bulk Ho helical pitch; we would like to 



emphasize that the allowance for at least six inter-layer interactions in the mode 
nian ([3]) turns out to be essential in order to be able to observe such behaviours 



Hamilto- 



B. Global Film Properties 

In this Section we will analyse some macroscopic thermodynamic quantities of the film, 
and for clarity reasons the attention will be again focused mainly on n = 12. We will show 
results pertinent to the magnetic specific heat, the chirality, Eq. (jHJ), and the average order 
parameter M denned in Eq. (jSJ). 
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FIG. 11: Specific heat C v vs. temperature for thickness n = 12 and lateral dimension 
L = 48 (star), 32 (triangle), 24 (diamond), 16 (square). 

The first quantity we consider is the specific heat. In Fig. [TT| C v at different L is displayed. 
Its behaviour clearly suggests the presence of two different phase transitions: In fact, two well 
separated maxima appear in Fig. [11] at T = 113.1 K and T = 119.4 K for L = 48, i.e. close 
to Tc i 2 (6) and Tc jl2 (l ... 5) respectively, making the maxima clear footprints of the block 



19 for n — 16 film: indeed, the 



phase regime. Such feature could not be observed in Ref. 
thinner temperature range where the block phase is present, joined with the broad character 
of the maxima for the finite-size samples investigated, made the two maxima coalesce and 
impossible to be resolved. Therefore, differently from what happens for n = 12, in that case 
the magnetic entropy seems completely released around Tc,i6(l ... 4, 7 ... 10, 13 .. . 16). 

The onset of a HM/Fan configuration along the perpendicular film direction can be probed 
by looking at a related quantity like the chirality, which is plotted in Fig. [12] together with its 
susceptibility. Xk(T) does not show any anomaly in the proximity of the highest temperature 
maximum of C v , i.e. around 120 K, while a clear peak appears at T ~ 114.2 K, becoming 
more and more sharp as L increases. In order to estimate the transition temperature Tjv(12), 
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FIG. 12: (a) Chirality vs. temperature at n = 12. The normalization factor is fixed by the 
bulk factor sinQ^ (see Eq. ([6])); error bars are smaller than point size. The continuum lines in 
the inset are obtained by multiple-histogram algorithm, (b) Xk obtained by multiple-histogram 
algorithm. The largest relative error is 0.5% for L = 48. 

a finite size scaling analysis of the quantities defined in the Eqs. (JH]), (|TU|) . and (ITT]) , with 
O = k, has been carried out. In Fig. [T3h <9s(k) is reported to show the typical behaviour of 
such quantities, while Fig. [T5b shows details of the fitting and extrapolation procedure. We 



obtain v using it as a free fit parameter in the equation 4fJ, |45|: 



In =L 1/u X (tL 1/u ) 



(14) 



where t = ^^7^, and X is an opportune scaling function. At the phase transition T/v(12) 
(i.e. for t = 0) we can consider the scaling relation ffT4l) as (dp In ^) max oc L x l v . Therefore, 
through an adequate fit, shown in the inset of Fig.[T3b. we have obtained v = 0.79(2). Using 
such value for z/, we can estimate T/v(12) from Eq. f[T2"|) by looking at the temperatures where 
(Xk)> and dp (In k) acquire their maximum value and extrapolating them against l}l v 

as shown in Fig. [T3"b . The final result obtained is T/v(12) = 113.6(1) K, a value definitely 
comparable with Tc,i2(6) = 113.4(4) K: we are thus lead to conclude that the onset of a 
helical/fan order in the film is only possible when all layers order, so that it is the last 
ordering spin layer, the 6 th one for n — 12, that drives the overall film transition to HM 
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FIG. 13: (a): dp{n) vs. temperature for L =12,16,24,32,40,48 (symbols as in Fig. EJ. (b): 
(12) plotted vs. L^ 1 ^ obtained by finite size scaling extrapolation of the three observables, 
with fitted v value. Inset: plot of the maximum value of In [dp (In k)\ vs. In L together with its best 
fit function (see Eq. ([P 



order. 

An issue largely debated in literature 
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44 



50 



5l| concerns the order of the chiral 



transition: In our MCS, in the whole thickness range here analyzed, we did not observe any 
double-peaked structure in the equilibrium energy distribution at Tjv(n), i.e., no explicit 
indication for a first order phase transition is given by our investigation. Anyway a first 
order transition can not be completely excluded, as suggested in Ref. |44j, where the author 



reasonably observes that a firm evidence for a first order transition can be obtained only 



44J 



when the sample is much larger than the largest correlation length 

The average order parameter M denned in Eq. (jSJ) turns out to probe the physical prop- 
erties of the system in a way more similar to what is done by the specific heat than by k, 
bearing it signatures of the onset of both spin block and HM/Fan phases, as it is apparent 
by looking at the related quantity xm- As an example, xm at n = 9 is reported in Fig. [T4l 
two anomalies are present at T ~ 112.0K and T ~ 117.0 K, i.e. at temperatures roughly 
corresponding to 7jy(9) = 111.2(5) K and Tc t g(l . . .4) = 115.9(4) K. We may observe that 
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FIG. 14: Susceptibility xm at n = 9 and lateral dimensions L = 24 (diamond), 32 (circle), 40 
(up-triangle) , 48 (square), and 60 (down-triangle) 

while the qualitative behaviour of xm is similar to that of the specific heat, the peaks in xm 
are sharper and display the finite size scaling typical of a critical quantity, thus making it a 
better probe to locate transition temperatures. 

As the film thickness decreases, the chirality and its related observables display a be- 
haviour similar to that at n = 12 up to n = 9, despite a shift from HM to fan-like order (see 
Fig. Hf); a largely different qualitative behaviour is instead obtained for n < 8, as already 
discussed in Sec. IV Al 

VI. DISCUSSION AND CONCLUSION 

In this paper the magnetic properties of thin Ho films have been carefully investigated by 
extensive MCS, assuming the model pertinent to bulk structure. Regarding the magnetic 
order below T/v(n), it has been showed as, by decreasing the number n of spin layers building 
the film, a progressive rearrangement of layer's magnetization from a helical to a fan-like 
structure (i.e. S XX {Q) ^ S yy (Q)), and finally to an essentially FM order for n < 7, is 
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observed. Moreover, for film thickness n = 8 the structure factor analysis has clearly revealed 
that once a finite magnetization has been established in every layer, an FM layer arrangement 
firstly appears which transforms to a fan-like configuration as the temperature is further 
reduced, as shown in Fig. [TU1 

Above all that, the system presents very interesting properties around the critical region 
when the film thickness is comparable with the bulk holmium helical pitch, i.e. for 9 < n < 
16. A spin block phase regime is observed in a wide range of intermediate temperatures. In 
this window some inner planes (n = 16 has a more complex block structure, as discussed 



in Ref. 



191 ) are in a paramagnetic configuration, while the other ones, close to the surfaces, 



appear to be in a quasi-FM ordered state (for example, when n = 12 we have obtained a 
spin block configuration where the magnetization rotate of an angle A<pi ~ 10° when moving 
from one spin layer to a neighbouring one within the same block, Fig. [7]). Everytime a spin 
block configuration appears, neighbouring ordered blocks line up in an antiferromagnetic 
way. What's more, also the study of macroscopic thermodynamic quantities, as the total 
energy, the order parameter defined in Eq. (JTj) and their derivatives, confirms the presence 
of such large critical regions. 

It is worth to remark that making use of all the six inter-layer coupling constants exper- 



imentally deduced by Bohr et al. in Ref. [20|, it is seen that the competition among surface 
effects and frustrated inter-layer interactions do not entail a simple adjustment of the surface 
planes only, but the magnetic critical properties of the whole film are strongly modified as 
well. Moreover, the results here presented, while confir ming that most of the predictions 
of the MFA employed by Jensen and Bennemann in Ref. Il3| are qualitatively correct, also 
show unambiguously that the thermal fluctuations play an essential role, so that the ability 
to include their effects is extremely important in order to have a full comprehension of the 
block phase phenomenology in these films. 

A detailed study of the chiraliy, Eq. ([6]), has shown that k correctly describes the estab- 
lishing of a global helical/fan order at T/v(n) for n > 9, but such quantity does not result 
critical in the temperature region where the spin block phase structure appears. Such be- 
haviour of k is also observed in the borderline case n = 8, where k does not present any 
anomaly at the single-layer ordering temperature Tc${\ ... 8). 

Another important issue of non-collinear film structures is the impossibility to describe 
T^f(n) vs. n through a well established scaling relation. In fact, as discussed in the Intro- 
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duction, making use of the empirical relation ([I]) one can easily only locate the thickness no 
which signalizes the disappearance of the HM order. From our results we can reasonably 



estimate no = 7 4- 9: The partial disagreement with the experimental results 15J (no ~ 10) 



could be a consequence both of defects and of the uncertainty with which the thickness of 



the Ho film experimental samples is known 241 ]. 



Finally, a possible comparison between the MCS results and the experimental data 15, 



161 ] would require particular care, in order to avoid naive considerations. Concerning the 



identification of the magnetic order type, we have shown [19J that the static structure factor 
alone is not able to distinguish between helical and spin block order. At the same time we 
think that the characterization of such spin block phases in non-collinear magnetic thin films 
could be very useful for experimental future works, being the coupling constants strongly 
dependent both on the used deposition substrate and on the employed deposition technique 
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521. 
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